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A stationary viscous jet falling from an oriented nozzle onto a moving surface is studied, 
both theoretically and experimentally. We distinguish three flow regimes and classify 
them by the convexity of the jet shape (concave, vertical and convex). The fluid is 
modeled as a Newtonian fluid, and the model for the flow includes viscous effects, inertia 
and gravity. By studying the characteristics of the conservation of momentum for a 
dynamic jet, the boundary conditions for each flow regime are derived, and the flow 
regimes are characterized in terms of the process and material parameters. The model is 
solved by a transformation into an algebraic equation. We make a comparison between 
the model and experiments, and obtain qualitative agreement. 



1. Introduction 

During the fall under gravity of a viscous jet from an oriented nozzle onto a moving 
surface, called the belt, three flow regimes can be distinguished. The flow regimes are 
characterized by the jet shape and depend on the process parameters such as dynamic 
viscosity of the Newtonian fluid, flow velocity at the nozzle, belt velocity, and falling 
height. 

In the first flow regime, apart from a boundary layer at the belt, the jet shape is 
concave and aligned with the nozzle orientation, and resembles a ballistic trajectory; see 
Figures 2(a) and |2(b)| We call this flow regime concave and the associated jet a concave 
jet. The concave jet occurs for large flow velocity at the nozzle and small viscosity. 

The jet in the second flow regime has a straight, vertical shape, apart from boundary 



layers at the nozzle and at the belt; see Figure 2(c) We call this flow regime vertical, and 



the associated jet a vertical jet. The vertical jet happens for large falling heights, large 
viscosity, and small flow velocity at the nozzle. 

In the third flow regime, apart from a boundary layer at the nozzle, the jet shape 
is convex, and the jet touches the belt tangentially; see Figure [2(d)] This flow regime 
we call convex, and the associated jet a convex jet. The convex jet occurs for high fluid 
viscosity, large belt velocity, small velocity at the nozzle, and small falling height. 

This paper is a continuation and a generalization of our previous work on the convex 



jet flfflod et q?.|p007| ). 

The fall of viscous jets or sheets from a nozzle oriented vertically down onto a fixed 
surface has been widely studied. Here one can observe unstable behavior; see |Taylor 



(1969 



Skorobogatiy fc Mahadevan (2000), Yarin & B.M. (19961, Ribe (20031, Ribe 



(20041, Cruickshank (19801 and Tchavdarov et al. (19931. Vertically falling viscous jets 
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have b een studied in |Clarke| ( |1966[ ) , |Clarke| ( |1968[ ) , |Adachi| ( |1987[ ) and |Sauter fc Buggischj 
(2005). Experimental investigations of steady and unsteady flows of a viscous jet falling 
under gravity onto a moving surface from a vertical nozzle were presented in Chiu-| 



Webster & Lister (20061 and Morris et al. (20081. In Chiu- Webster & Lister (20061 and 



Ribe et al. ( 2006 1 the steady flow is modeled and the parameter region of the steady flow 



is determined in terms of the falling height and the surface velocity. 

However, the previous publications make no distinction between concave and vertical 
flows, and because the nozzle is oriented vertically down the concave flow is not recog- 
nized as a separate regime. In this paper, we fully describe all three flow regimes. To 
describe the jet we use a model which includes effects of inertia, viscosity and gravity, 
but neglects surface tension, bending stiffness and air drag. The fluid is considered to be 
incompressible, Newtonian, and temperature effects are neglected. We allow the nozzle 
orientation to vary between horizontal and vertically down. By studying the characteris- 
tics of the equation of momentum conservation, we determine the parameter regions for 
each flow regime. Consideration of the characteristics as being the directions of informa- 
tion propagation explains why and when each of the three flow regimes occurs and gives 
the correct boundary conditions for each flow regime. To validate our theoretical results 
we perform experiments of the jet falling from the oriented nozzle onto the moving belt. 
We find a qualitative agreement between the experimentally observed and the theoretical 
values of the positions of the touchdown points for different belt velocities. The model 
presented in this paper can also be used to describe the fall of viscous sheets onto a 
moving surface. 

The structure of the paper is as follows: In Section [2j we describe the experiments 
of the fall of the viscous jet onto a moving belt, and present the experimental results. 
In Section [3] the model equations are derived and simplified to a first-order differential 
equation on unknown domain. The analysis of the characteristics of the conservation of 
momentum equation for dynamic jets in order to derive correct boundary conditions is 
given in Section]?] In Section [5] we present some results from the model, and in Section [6] 
we compare them with experiments. The characteristic features of the three flow regimes 
are summarized in[7j and some conclusions are made in Section [8] 



2. Experiments 

In this section we describe experiments of the fall of a thin jet of a Newtonian fluid 
onto a moving belt. We focus on the shape of the jet between the nozzle and the belt. We 
describe the experimental setup, report our observations and present some conclusions 
from the experiments. 

2.1. Experimental method 

A viscous fluid, polybutene Indopol H-100, is pumped to a nozzle and allowed to fall from 
the nozzle onto a moving belt; see Figure [T] The belt is wrapped around two horizontal 
cylinders at the same height. The left cylinder is connected to an electric motor, to move 
the horizontal belt from the left to the right with a constant speed. 

The nozzle is placed above the belt. The nozzle - belt distance and the belt and 
the nozzle orientation can be varied. A screw pump producing a constant flow rate is 
connected to the nozzle. The flow rate was measured by weighing the fluid collected 
from the nozzle during 30 s. In the experiments, two different nozzles were used, with 
diameters of 1mm and 0.4mm. 

The experimental setup allows us to change the nozzle position and orientation, belt 
velocity, and flow velocity from the nozzle. For all experiments the same fluid is used. 
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Parameter name 





Value 


Unit 


"belt 


- 5 


m/s 


V nozzle 


0.4 - 1.2 


m/s 


L 


0.01 - 0.07 


m 


Q^nozzlc 


-9 - 38° 




V 


0.047 


m 2 /s 


P 


880 


kg/m 3 


^nozzle 


1 or 0.4 


mm 



belt velocity 

flow velocity at nozzle 

distance between belt and nozzle 

nozzle orientatiorjf] 

kinematic viscosity of fluid 

fluid density 

nozzle diameter 



t The angle between the nozzle orientation and the 
horizontal direction, positive for downwards-pointing 
nozzle. 

Table 1. Values of the experimental parameters 




Figure 1. The scheme of the experimental setup for a jet falling onto a moving belt. 



The fluid is considered to be Newtonian. No nonlinear effects such as die swell near the 
nozzle were observed. The values of the experimental parameters are given in Table [T] 

2.2. Experimental results 

First, we describe a typical sequence of experiments. We start with Vbeit close to zero 
and make sure that L and v nozz \ e are chosen such that the shape of the jet is concave, 
resembling a ballistic trajectory; see Figure 2(a) To obtain the concave jet shape the 
nozzle should not point down vertically, and therefore we put a nozz \ c < ir/2. Next, we 
gradually increase Vbeit and study the evolution of the jet shape. 

For small Wboit the jet shape is concave with an unstationary region near the belt; 
see Figure [2(a)| By increasing Ubcit, we observe that the unstationary region near the 
belt transforms into a stable bending region where the jet bends to the horizontal belt 
direction; see Figure |2(b)| The jet shape in this region resembles the backward-pointing 
heel, reported for the vertically falling jet in Chiu- Webster & Lister (20061. 



When we increase Vbeit further, the jet shape approaches the vertical direction. Note 
that this direction of change is often considered counterintuitive. In this case, the contact 
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(a) Ubeit = 0.093 m/s 



(b) Vbeit = 0.750 m/s 




(c) Ubeit = 1.708 m/s 



(d) w b ei t = 3.200 m/s 



Figure 2. Jet shapes for the set of experiments with d nozz i c = 1mm, u nozz i c = 1.061 m/s, 
L = 0.054 m and a nozz i c — 37. 3°. T he belt m oves from the left to th e right. The jet shape 
changes from concav e (Figures 2(a) and 2(b) I to vertical (Figure 2(c) I, and from vertical to 
convex (Figure 2(d) I, while Wb e lt is increased. 



point with the belt approaches the vertical projection of the nozzle position. As a result, 
for fbcit large enough, the main part of the jet between the belt and the nozzle is purely 
vertical; see Figure 2(c) The bending region near the belt remains, and a new bending 
region near the nozzle appears. Near the nozzle the jet bends from the nozzle orientation 
to the vertical direction. 

Further increase in Vbcit results in the disappearing of the bending region near the belt. 
The jet shape becomes convex everywhere, except for a bending region near the nozzle; 
see Figure 2(d) The touchdown point moves away from the nozzle in the direction of the 
belt motion as Vbeit increases. 

Summarizing the results of the experiments, we observe a concave jet shape for small 
Ubeitj except for a small bending or unstable region near the belt. With increasing Ubeit 
the jet shape becomes vertical, except for small bending regions near the nozzle and the 
belt. Further increase of Vbeit leads to a convex jet shape, except for a small bending 
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Figure 3. Experimentally obtained positions of the touchdown point x en d for different Vhelt 
indicated by •. The other parameters are dnozzie = 0.4mm, v uozz \ e = 1.147m/s, L = 4.1cm and 

r o 

Q^nozzlc — O . 




region near the nozzle. This gives a characterization of the jet flow by its shape, i.e. 
concave, vertical and convex. 

A convenient way to quantitatively study the jet is to look at the horizontal position 
Xcnd (see Figure [3]) of the touchdown point at the belt, relative to the nozzle position. For 
the concave jet with Wboit small the contact point is away from the nozzle in the direction 
of the nozzle. With increasing Wbeit, the touchdown point moves first towards the nozzle 
position (x cn d decreases) until the jet becomes vertical (a; en d = 0), stays vertical for some 
time and then moves away from the nozzle position in the direction of the belt motion 
(a; on d increases). Figure [3] suggest that the first five dots are in the concave flow regime, 
the sixth, with x cn d — 0, represents a vertical flow, and the remaining ones (7th and 
higher) are in the convex flow regime. 

For the concave jet in Figures |2(a) and |2(b)[ where the nozzle diameter is 1 mm, the 
jet near the nozzle is aligned with the nozzle orientation. When the nozzle diameter is 
0.4 mm, the jet is not aligned with the nozzle orientation. This can be seen in Figure [5j 
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Figure 5. Experiment with d nozz i c = 0.4mm, u nozz i c = 1.188m/s, «belt = O.lm/s, L = 0.068m 
and a nozz io = 5°. The jet has a concave shape, which is not aligned with the nozzle orientation. 

In the next section we present a model that characterizes the flow type and, among 
other things, predicts the relation between a; on d and fbeit- 

3. Modeling and analysis 

In this section we present our model of the fall of a viscous jet onto a moving belt. To 
model the flow we use a thin-jet approximation and include effects of inertia, viscous ten- 
sion and gravity. We assume the fluid to be incompressible, isothermal, and Newtonian. 
We neglect surface tension, bending stiffness, and air drag. Therefore possible bending or 
buckling regions at the nozzle or at the belt are disregarded. The jet is described by the 
equations of conservation of mass and momentum. First we formulate the equations for 
the dynamic jet, which are used in Section [4]to justify or choice of boundary conditions. 
Next, we partly solve the stationary jet equations and make an analysis showing that the 
jet can have only three possible shapes: concave, vertical and convex. Finally, we refor- 
mulate the problem by deriving an equivalent algebraic equation which is convenient for 
further analysis. 

The jet is modeled as a curve in the x, z-plane of unknown length s on d,; see Figure [4j 
The curve is parameterized by its arc length s, with the origin s = at the nozzle and 
s — Send at the touchdown point at the belt. The position of a certain point s of the thin 
jet at time t is described by its position vector r = r(s, t) with respect to the origin 0, 
which is chosen at the nozzle point. 

A local coordinate system in a point s having as basis the tangent and normal vectors 
e t , e n , is constructed at each point of the jet. The angle between the tangent vector and 
horizontal direction is 0. The horizontal distance between the nozzle and the touchdown 
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point at the belt is x en d- The flow velocity in a point s of the jet is v — v(s, t). The jet 
at the touchdown point has the same velocity as the belt Wbcitj & n d the flow velocity at 
the nozzle is u noZ zie- 

The system of equations describing a thin dynamical jet in two dimensions can be 



found in variety of publications i.e. Roos et al. (20061, Yarin (19931, and lEntonov & 



Yarin (19801. It consists of the laws of conservation of mass and momentum 

At + (Av) s = 0, (3.1) 

pA{v tt + r s (v t + w s ) + v 2 r ss + 2vr st ) = P s r s + Pr ss + K', (3.2) 

respectively, where A = A(s, t) is the cross-sectional area, P = P(s, t) is the longitudinal 
force, and K' = K'(s, t) is the external force per unit of length of the jet. The longitudinal 
force P is given by a constitutive law, and in the case of a Newtonian viscous fluid it is 
equal to 

P = 3vpv s A. 

Finally for K' we take 

K' = pAg, 

the gravity force per unit of length of the jet (external air drag is neglected). 

The stationary versions of the equations ( |3.2| -(3.1 1 together with (3.3)-( [3~4| 
condition for s as the arc length are 

A(r s vv s + v 2 r ss ) = 3is(v s Ar s ) s + Ag, 

(Av) s = 0, 
IrJ = 1. 



(3.3) 

(3.4) 
and the 



(3.5) 
(3.6) 
(3.7) 



Thus, we have three differential equations, (3.5l-(3.7l, for the unknowns r, v and A. Next 
we describe the boundary conditions. 

For the velocity v we prescribe two boundary conditions: at s — 0, the flow velocity at 
the nozzle is 



v(0) = IWzlc, 

while at s = s cn( j the jet sticks to the belt, so 

f(Send) = ^belt- 



(3. 



(3.9) 



The boundary condition for A follows form the known cross-sectional area of the nozzle 



as 



AO) = T<& 



(3.10) 



The fixed vertical distance between the nozzle and the belt gives the additional constraint 

sin6(s)ds = L. (3.11) 



To make the system ( 3.5 )-( 3.11 ) complete we need two boundary conditions for r. Since 
the position r is with respect to the fixed nozzle, we have one boundary condition for r 

r(0,i) = O. (3.12) 

The second boundary condition is chosen later in this section and our choice is justified 
in Section |4j By integrating (3.6 1, using (3.8) and (3.101, we find that 

F 
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where the mass flux is given by F — pv nozz \ c nd^ ozzlc / A. We eliminate A from (3.51 to 
obtain 



r s v s + vr ss = 3v(r s v s /v) s + g/v. 
Next we introduce a new variable £ by 

v s 



•v 



(3.13) 



(3.14) 



which stands for the scaled momentum transfer through a jet cross-section and plays a 

(3.15) 



crucial role in our further analysis. By use of £ we write (3.131 as 



Using e t = r s , and (e t ) s — — s e n , we can write (3.15) in components as 

gsin(e) 



and 



e. 



gcos(6) 



(3.16) 



(3.17) 



Equation (3.17) requires a boundary condition for 0; this is related to the question of 
boundary conditions for r. 

We scale the system as follows: the length s is scaled with respect to 3^/i> nozz i e , and 
the velocity v with respect to iWzic- Then, (3.161, (3.171, (3.141, (3.8|, ( 3.9 > and (3.111 
become 



Asin(e) 

) 

V 

Acos(e) 



£ = v , 

v 

v(0)=l, 
v(s end ) = Dr, 



(3.18) 
(3.19) 

(3.20) 

(3.21) 
(3.22) 

sin(6(s))ds = Re. (3.23) 

Rc = v nozz i e L/(3v) is the Reynolds number, Dr = w bc it / Vozzie is 
the draw ratio, and the scaled s en d becomes s e ndi'nozzie/(3i / ). The dimcnsionlcss number 
A is related to the Froude number Fr = v nozz i c /\/gL and Re as A — l/(ReFr 2 ). After 
scaling the system is described in terms of three positive dimensionless numbers, which 
define a parameter space & as 

& = {(A,Re,Dr) : A > 0, Re > 0, Dr > 0}. (3.24) 

The nozzle orientation a nozz i e only appears in the boundary condition for for the 
concave jet (3.38), and is considered to be fixed. 



Here, A = 3gv/vl ozzlc , 



By replacing the material coordinate s by the time variable r, according to 

ds = v{r)dT, 



the system (3.18 1-( 3.23 1 becomes 

£ T = Asin(9), 



(3.25) 
(3.26) 
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(3.27) 

(3.28) 

(3.29) 
(3.30) 

(3.31) 



Here, T en d is the result of the coordinate transformation (3.251 of s cnc i = / °" d v(t)cLt. 

(3.32) 



Next, we solve (3.261 and (3.27), using the first integral 

£sin(6) = At+Ci, 
to obtain 



£ = ±V / A 2 t 2 + 2Ac 1 t + c 2 , 

n x- ■ ( Ar + a 
B = ± arcsm 



VA 2 



(3.33) 
(3.34) 



-2A Ci t + c 2/ 

Here, C\ and c 2 are unknown constants to be determined later. 

In the analysis we restrict ourselves to solutions with £ [0,7r/2]. Then, we conclude 
from ( 3.26 1 that £ is a strictly increasing function. Therefore, we distinguish three possible 
situations for the sign of £: always positive, a sign change from negative to positive, and 
always negative, i.e. 



< £(Q) < £(r cnd ), 
m < < £(r cnd ), 
m < £(Tcnd) < 0. 



(3.35) 
(3.36) 
(3.37) 



If (3.351 holds, then it follows from ( 3.27[ ) that G is a strictly increasing function for 
O < 7r/2, implying that the jet has a concave shape. As will be justified in Section [4] we 
prescribe the nozzle orientation angle as the boundary condition for 0, i.e. 



9(0) 



Substitution of (3.381 into (3.33 1- (3. 34 1 gives 



e 



A 2 t 2 + 2A^/c 2 sin(a nozzlc )r + c 2 , 



At + v /c^sin(a nozz i ) 



y/ A 2 t 2 + 2A A /c^sin(a nozz i )r + c 2 



(3.38) 

(3.39) 
(3.40) 



Because (3.351 implies a concave shape, we refer to a jet satisfying (3.351 as a convex jet. 
For (3.36) to hold, there must exist a t* G [0,T on( j] such that £(r*) = 0. Then from 

(3.41) 



fl3.27| , it follows that 6(r*) = tt/2. Substituting r* into ( |3.34[ ), we have 

At 



Cl 



(T~)* + 2AC\T* + c 2 

giving = C2. This implies that 

6 = tt/2, 

for all t G [0,T on( j], and hence the jet is vertical, and 

£ = At + Ci . 



1. 



(3.42) 
(3.43) 
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For £ obeying (3.361, we obtain £(t) = At — \J (c-i) 2 ■ Because (3.36) implies a vertical 
shape, we refer to a jet satisfying (3.361 as a vertical jet. Note that for the vertical jet, 



as will be shown in Section [4] no boundary condition for is necessary. 

If ( |3.37 1 holds, then it follows from (3.27) that is a strictly decreasing function for 
< 7r/2. In this case the jet has a convex shape. As will be justified in Section [4] we 
require tangency for the jet at the belt, i.e. 



Then 



0(T c „d) = 0. 

i = - V / ^ 2 T(r-2r cnd ) + c 2 , 

4(T cnd - T) 

^A 2 t(t - 2r end ) + c 2 



= arcsin 



(3.44) 

(3.45) 
(3.46) 



Because (3.371 implies a convex shape we call a jet for which (3.371 holds a convex jet. 



By substituting the found solutions for £ and into ( 3.28 1-( 3.31 1 for the three situa- 



tions (3.35 !-( 3.37) we successively obtain 



VA 2 r(r-2r cnd )/ U ; 2 + l 



v(Q) 

v(T cnd ) 



Re = 




2Arwsin(a nozz i e ) concave jet, 
vertical jet, 
convex jet, 



At-\-w sin(a n 



e)_ 



'° -J A 2 t 2 +w 2 +2Atw sin(a, 

f^v(r)dr 



= v (t ) dr concave j et , 

vertical jet, 
convex jet, 



(3.47) 

(3.48) 
(3.49) 

(3.50) 



y/ A 2 T(T-2T Bnd )+W 2 

where w = £(0). We refer to the situations of concave, vertical and convex jets as concave, 
vertical and convex flow regimes, respectively. 



For given w G R and flow regime, the problem ( 3.47 1-( 3.49 1 has a solution v(t;w) and 
T e nd(w), where T cn d(w) satisfies (3.491. Here, we assume that for any w, ( 3.49[ ) has only 
one solution, which is not always true. However, this allows us to illustrate a solution 
procedure. 



Substituting v(t; w) and T en< x(w) into the integrals (3.501, we obtain the functions of w: 
rTondW 4r + wsin(a nozz io) 



^conc(w) 



icomr(w) 



o ^fA 

T en d(l") 



II 



t 2 + w 2 + 2Atw sin(a noZ zic 
Wvort(T; w)dr 

A(T end (w) - t) 



y/A 2 T(T-2T end ( W ))+W 2 



v CO nc(T~; w)dr concave jet, 
vertical jet, 
Vconv (t ; w)dr convex jet. 



(3.51) 

Here, we denote by v COIic (T\ w), v vcrt (T;w) and v conv (T;w) the solution of (3.471 for a 
concave, vertical and convex jet, respectively. According to ( 3.35 >-( 3.37 > , I C onv(w) and 
Ivert(w) are defined for w ^ 0, and I conc (w) for w > 0. With (3.51 1, solving (3.47l-( [3~50 l 
is equivalent to solving the algebraic equation 

J 7 (w) = Re, (3.52) 

where ? stands for an unknown jet flow regime. Therefore, a study of existence and 
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Figure 6. Parameter regions for three flow regimes ^ conc , & ve rt and ^ c(mT . The border 
between ^ conc and ^ VO rt is Si and the border between ^ conv and & VCI t is 52- 



uniqueness of a jet solution results into a study of the existence and uniqueness of a 
solution to the algebraic equation (3.521. 

At this point, we like to briefly recapitulate the main steps in our solution procedure. 
We do this, as an example for the concave flow; the other cases are completely analogous. 
The steps are: 

(a) Solve v = v conc (T; w) from ( 3.47 1 x , with use of the boundary condition (3.481. 
(6) Find r ond (w) from ( |3.49| as Vc onc(T cnd (w); w) = Dr. 

(c) Calculate I CO nc(w) from (3.511. 

(d) Solve w from < pT52 |. 

The partitioning of the parameter space into the regions of concave ^ C onc, vertical 
^vort and convex # conv jets is presented in Figure [6| The partitioning follows from the 
solutions of (3.47 )-( 3.50 ) with the additional condition £ (r en d) = for the border between 
invert and ^conv, and £(0) =0 for the border between ^ conc and ^vert- 

For any set of parameters from <^ conc , or ^ V ort , or # conv , a corresponding solution of 
(3.47 1-( 3.50 1 exists. This solution is unique in the case of vertical and convex flow. For 
the concave flow a solution might not be unique when the nozzle does not point vertically 
down and Dr > 1. In this case up to two solutions exist for a concave jet and one solution 
for a convex or vertical jet. The questions of existence and uniqueness are fully treated in 
the upcoming paper [Hlod fc Peletier| ( p009[ ) , and for the convex jet in |Hlod et «L| ^2007). 

In the next section we justify our choice of the boundary conditions for 9, i.e. (3.38l 
and (3.441, for concave and convex flow, and we explain that no boundary condition is 
needed for 6 in case of vertical flow. 



4. Justification of boundary conditions for 

In this section we explain our choice of boundary conditions for made in the previous 
section. We use the highest-order part, or principal part, of the conservation of momentum 
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equation (3.2 I 



r ti + 2vr st + v£,r ss = [r tt + +2vr s 



v 2 r„ 



v.r. 



(4.1) 



This equation is of hyperbolic type in the neighborhood of s = s eru j if £(s C nd) is close to 
zero, and in the neighborhood of s = 0, if £(0) is close to zero, because v s — v2 + £v > 
in theses situations. We consider the direction of the characteristics of (4.1), either at 
s = s on d if £(s C nd) changes sign, or at s — if £(0) changes sign. This directly yields 



the number of boundary conditions, which must be prescribed at s 



„d or s = 0. The 



reason is that the number of boundary conditions at any point of the boundary is equal 
to the number of characteristics pointing inside the domain at this point; (see |Godlewski| 



fc Raviart||1996 p 417) 



The characteristics equation (see Davis 2000 p 57) for (4.1 1 is 



2vz- 



(4.2) 



where z represents the velocity of a characteristic. Equation (4.2) has the solutions 



~2 



(4.3) 



According (4.3| and (3.20), the directions of the characteristics of (4.1l depend on the 



sign of £ as follows: 

(a) If £ > then z\ > and z-i > 0, i.e. the two characteristics point to the right. 

(b) If £ = then Z\ > and z^ — 0, i.e. one characteristic points to the right and one 
is stationary. 

(c) If £ < then z\ > and Z2 < 0, i.e. one characteristic points to the left and one 
to the right. 



From the characterization of the flow regimes ( 3.35 )-( 3.37) we infer that: 

• At s — two boundary conditions for r(s, t) are necessary in case of a concave jet 
(£(0) > 0), and only one in case of a vertical or convex jet (£(0) < 0). 

• In case of a convex jet, one boundary condition for r is necessary at s = s C nd 
(C(send) < 0), and no in case of vertical or concave jets (£(s en d) ^ 0). 

For all three situations we prescribe the nozzle position. In addition, for the concave jet 
the nozzle orientation is prescribed by (3.381, and for the convex jet the tangency con- 



dition (3.44 



(3.381 and 



with the belt is prescribed. This justifies our choice of boundary conditions 



3.44) for the stationary problem. 



The analysis of characteristics, as directions of information propagation, explains why 
the nozzle orientation influences the jet shape only in the case of concave flow, and why 
the belt orientation influences the jet shape only in case of convex flow. 

• In concave flow all information about the jet shape travels from the nozzle to the 
belt. Therefore, not only nozzle position but also nozzle orientation is relevant for the 
jet. In addition, no information on angle travels back from the belt. 

• In vertical flow only one characteristic (at the nozzle) points inside the domain. 
Therefore, no information about nozzle orientation or belt movement direction influences 
the jet shape. Thus, in vertical flow the nozzle and the belt orientations are irrelevant 
for the jet. 

• In convex flow one characteristic points inside the domain at the nozzle and one at 
the belt. Hence, information about the direction of the belt movement influences the jet 
shape, and therefore, the belt orientation becomes relevant in convex flow. 
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All flow regimes are possible 



1. ^4* 3. 4. 5. 6. 
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Figure 7. Regions of A and Re for which, for all Dr, different flow regime exist. In the region 
{A < A*, Re > ,'%i(A)}, only concave jets are possible, in {A > A* , Re > .^(A)} only vertical 
jets, while in {A > 0, Re < ^s(A)} all flow regimes are possible. In the remaining region 
between Mi and concave and vertical jets, but no convex jets, are possible. 

5. Results from the model 

In this section we present some results from our model. We analyse the partitioning of 
the parameter space. Next, we investigate changes of the flow type if one of the physical 
parameters (L, v, Ubeitj ^nozzle) is varied. We describe the trajectories of the process 
parameters in the parameter space & , and we illustrate the jet shape evolution. Note 
that the only possible transitions between flow types are between ^vert and # COI i V , and 
between ^Vert and ^ conc ; see Figure [6] 

The projection of the regions for the three flow regimes onto the (A, Re)-plane is 
depicted in Figure [7] and following is valid for all values of Dr. We observe a region {A < 
A*, Re > Mi(A)} where the jet is concave, and a region {A > A*, Re > Mi(A)} where 
it is vertical. In the region between £%2 and 8%§ concave or vertical flow is possible, 
but there can be no convex flow. Finally, in the region {A > 0, Re < ^(A)}, all three 
flow regimes are possible. Hence, it is only in the latter region, where Re < ^^{A) < 1, 
that a convex jet can occur. 

The parameter regions projection onto the (A : Dr)-plane is depicted in Figure [8] and 
following holds for all Re. We observe that for {A > A* , Dr < 1} only a concave jet is 
possible, while concave or vertical flow is possible for {A > A*, Dr < 1}. In the region 
{A > A* , Dr > @(A)} convex or vertical flow is possible, while in the rest of region 
{Dr > 1} all three flow regimes are possible. Hence, a convex jet can only occur if 
Dr > 1. 

Next, we study the evolution of the jet if one of the physical parameters varies as to 
change the flow type from convex to vertical. For a reference configuration we take the 
physical parameters L — 1cm, v = 0.047 m /s, Wbeit = 1.4 m/s, and u n0 zzio = lm/s, for 
which the jet is convex. Then, if we increase L, decrease v, decrease fbdt, or increase 
^nozzle, eventually the jet flow changes from convex to vertical. The corresponding curves 
in the parameter space £P are indicated in Figure [9] 
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All flow regimes 
are possible 


Vertical or convex jet 




Only concave jet 


Concave or vertical jet 



u.u — i 

0. I. A 3. 4. 5. 6. 

A 

Figure 8. Regions of A and Dr for which for all Re different flow regimes exist. In the region 
{A < A* , Dr < 1} only concave jets are possible, in {A > A* , Dr < 1} concave and vertical jets 
(no convex jets), and in {A > A* , Dr > S>(A)} convex and vertical jets (no concave jets). In 
the remaining region between Dr = 1 and @, all flow regimes are possible. 




Figure 9. Traces of the point (A, Re, Dr) as we change one of the physical parameters (L 
increases, v decreases, Vbelt decreases, u n ozzic increases). The curves originate at the same point 
in the region # collv and eventually leave ^ conv by crossing the separating gray surface at the 
points indicated by the dots. 
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(a) Convex jet shapes for different L: 1, 1.7, (b) Convex jet shapes for different v: 0.047, 
2.2, and 3.5 cm. The shape approaches the 0.026, 0.015, and 0.012 m 2 /s. The shape ap- 
vertical as L increases. proaches the vertical as v decreases. 



v (cm) 




x (cm) 




x (cm) 



(c) Convex jet shapes for different Ubcit: 1-4, (d) Convex jet shapes for different v noz2 ,i a : 
1.21, 1.11, and 1.08 m/s. The shape ap- 1, 1.16, 1.24, and 1.26 m/s. The shape ap- 
proaches the vertical as Ubcit decreases. proaches the vertical as ^nozzle increases. 

Figure 10. Shapes of the convex jet for different values of L, v, Wboit and fnozzic- The reference 
values are L = 1cm, v — 0.047 m 2 /s, fbcit = 1.4m/s, and « n0 z Z i c = lm/s. 



Changes of the jet shape while only one of the physical parameters L, v, Ubeit) or 
^nozzle varies as described above are shown in Figures 10(a)| |10(b)| |10(c)[ and |10(d)[ 
respectively. In Figures [9] and [10] we see that if the point (A, Re, Dr) approaches the 
boundary of # COIlv , the jet shape becomes vertical. If (^4,Re,Dr) is very close to the 
boundary of ^ conv the jet shape is almost vertical, except for the small region near the 
belt where the jet rapidly bends to the horizontal belt direction. 

The analysis of the parameter region for the concave jet is more complex than that 
for the convex jet. In case {A > A*, Dr < 1} ("Only concave jet" in Figure [8|, the 
flow is concave for all L. Similar, if {A < A* , Re > £%i(A)} ("Only concave jet" in 
Figure [7| , the flow is concave for all Vbcit ■ In a situation when v decreases or i> n0 zzic 



increases, A approaches zero and Re approaches infinity since A — 3gv/v^ 



zlc ' 



and 



Re = v 'nozzle L I (3f) ■ Thus, eventually the point (A, Re) enters the "Only concave jet" 
region in Figure [7] Hence, if the jet is not in the concave flow regime, decreasing v or 
increasing w ri ozzic makes the jet to become concave eventually. 

To illustrate the change of flow from concave to vertical, while only one of the param- 
eters L, i>beit, and w n0 zzic varies, we take the reference values L — 30cm, v = 0.2m 2 /s, 
^beit = 2 m/s, and i> n0 zzie = 1.5 m/s. Then, if we decrease L, increase v, increase Ubcit, or 
decrease ^nozzle eventually the jet flow changes from concave to vertical. The curves in 
the parameter space are indicated in Figure [TT| 

Changes of the jet shape for a noZ zic = tt/4, while only one of the physical parameters L, 
v, Wbeltj or ^nozzle varies as described above are shown in Figures [l^^p7T^[l2(cY| and 
|12(d)| respectively. In Figures 11 and 12 we see that if the point (A, Re,Dr) approaches 



1G 




Figure 1 1 . Curves in the parameter space as we change one of the parameters (L decreases, v 
increases, «beit increases, ii n ozzio decreases). The curves originates at the same point in the region 
J^conc and eventually leave ^ conc by crossing the separating surface at the points indicated by 
the dots. 

the boundary of ^ conc , the jet shape becomes more vertical. If (A,Re,Dr) is very close 
to the boundary of ^ C onc the jet shape is almost vertical except for the small region near 
the belt where the jet rapidly bends from the nozzle direction to an almost vertical one. 

The analysis above shows that the transition between the convex and the concave flow 
regimes as parameters continuously vary is only possible via the vertical flow. In the next 
section we compare the results from our model with those from our experiments. 



6. Comparison between the model and experiments 

In this section we validate our model using the results of the experiments described in 
Section [2] by comparing the corresponding relations between z ond and "belt- We compare 
the shapes from the experiments and the model, and discuss differences and similarities 
for jets in convex and concave flow regimes. 

6.1. Comparison of x en d 
We start by comparing the model predictions of the touchdown point with the ones 
obtained from the experiments in Figure [3] The following values of parameters are used 
"nozzle = 1.147 m/s, L = 4.1cm, v = 0.047 m 2 /s, d nozzlc = 0.4mm, and a nozz i c = 5°. 
In Figure [13] we present the relations between x cn d and "belt obtained from the model 
(solid curve) and the experiments (dots). The partitioning of the parameter space 3? 
gives the regions of "belt for the three flow regimes as indicated by the vertical lines (at 
"belt = "l and V2). For small i>bcit the flow is concave, for increasing t>bcit vertical, and 
for even larger Wbeit convex. 

A curve C(a; en d, "belt), describing the relation between a; en d and "belt as obtained from 



the model (solid curve in Figure 13 1 consists of three parts: 
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(a) Concave jet shapes for different L: 30, (b) Concave jet shapes for different v: 0.2, 
18, 13, and 12 cm. The shape approaches 0.26, 0.3, and 0.32 m 2 /s. The shape ap- 
the vertical as L decreases. proaches the vertical as v increases. 



y (cm) V ( cm ) 




(c) Concave jet shapes for different «belt : (d) Concave jet shapes for different v nom i e : 
2, 2.57, 2.86, and 2.95 m/s. The shape ap- 1.5, 1.4, 1.36, and 1.34 m/s. The shape ap- 
proaches the vertical as Ubdt increases proaches the vertical as w noZ zie decreases. 



Figure 12. Shapes of the concave jet for different values of L, v, Wbcit, ^nozzle- The reference 
values are L = 30 cm, v — 0.2m 2 /s, Wboit = 2 m/s, and Wnozzle = 1.5 m/s. The nozzle orientation 

IS Qnozzlc = 7I"/4. 



Concave jet part. For Uheit close to 0, x en ^ increases as Vheit increases, till C(x cn d, fboit) 
reaches its maximum. For further increasing Wbeit, a^ond decreases till C(a; e nd, ^bcit) be- 
comes vertical (at Vbcit = ^2)- After this point C(x on d, ^b it) bends back and x en d decreases 
to zero; during this phase i>beit decreases to v\. At this point where a; cn d becomes zero, 
the jet becomes vertical; in the preceding part of C(x en d, ^beit) the jet is concave. Hence, 
for t>boit € [Oj^] and a: e nd > 0, the jet is concave. 

Vertical jet part. This part of C(x en d) ^beit) is horizontal with x on d = 0, while i>bdt 
increases from v\ to t>2- In this part, the flow is vertical. 

Convex jet part. This part of C(a; Gn d, i^cit) starts at a; on d = and Vbeit = ^2, after 
which both a: C nd and Wboit increase. In this part the jet is convex. 

Looking at the shape C(x en d, ^bcit) in Figure 13 within the region v\ < Vhcit < «2, we 
notice that the solution is non-unique there: two concave and one vertical solution exist 
there. This illustrates the non-uniqueness of the jet solution for a nozz \ c < tt/2. 
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Xend/L 




Vbeit (m/s) 



Figure 13. Comparison of the relations between Ubcit and a; on d as obtained from the model 
C(x en d, Vbeit) (solid line) and from the experiment (dots) for L = 4.1cm, v = 0.047 m 2 /s, 
^nozzle = 1.147m/s, d n ozzic = 0.4mm, and a n ozzic = 5°. A vertical jet occurs for «b e lt G [«i,«2], 
and x cnd = 0. 



Figure 14. The regions of Vbeit and L for the three flow regimes and the experimental values of 
Zend- Here, two values of L are used 4.1 cm and 6.8cm. The other parameters are v — 0.047 m 2 /s, 
^nozzle = 1.147 m/s, dazzle = 0.4mm, and a n ozzic = 5°. The value of x cnd is small for u bc it in 
the parameter region for the vertical jet and increases if i>b e it goes away from this region. 

The experimental results have qualitatively the same tendency as the theoretical ones 
from our model. For Ubeit close to zero (the first three experimental values) the small 
maximum, or a kind of plateau, is found too and after that the values monotonically 
decrease till a; en d becomes zero. This represents the transition from the concave to the 
vertical jet regime. The following observation points all lie in the convex jet regime and 
they show a monotonic increase of x on( j with Vhcit- Hence, the behavior of the experimental 
data agrees, in general, with that predicted by the model. The only difference is (non- 
)monotonicity in the small region between v\ and v^. 

Although the theoretical and experimental results agree in a qualitative sense, quan- 
titatively significant differences are found. The values of x onc i predicted by the model 
for convex and concave flows are larger than the values obtained experimentally We 
comment on this in the next section 

Finally, we compare the predictions of the parameter regions in the (i>b e it) £)-plane for 




L (cm) 



4 
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2 







Three viscous jet flow regimes 19 




Figure 15. Comparison of the concave the jet shapes obtained theoretically and 
experimentally for L — 4.1cm, v = 0.047 m 2 /s, u n ozzie = 1.147m/s, «belt = 0.81 m/s, and 

ro 

^nozzle " • 



V (cm) 




Figure 16. Comparison of the convex jet shapes obtained theoretically and experimentally for 
L — 4.1cm, v = 0.047 m 2 /s, fnozzic = 1.147m/s, Vhelt = 4.398m/s, and a n ozzie = 5°. 

flow regimes obtained from the model with those from the experiments. To determine the 
type of flow regime from the experimental data, we observe the behavior of x cn( j as Hbeit 
increases: for concave flow x cn a decreases, for vertical flow a; en( j is close to zero, and for 
convex flow x on d increases. The results for two values of L are presented in Figure [l4j For 
the smaller L we clearly see that in the concave jet region x on( j decreases, stays close to 
zero in the vertical jet region, and increases in the convex jet region as i>bcit increases. A 
similar behavior of a; cn d is observed for the larger L, except in the convex jet region where 
the increase of x en d is less significant. Summarizing, we conclude that the experimental 
results for the parameter regions for the three flow regimes in the (ubeltj £)-plane agree 
with the theoretical ones. 

6.2. Comparison of jet shapes 

Theoretically and experimentally obtained jet shapes are presented in Figures [T5| and [T6] 
for a concave and a convex jet, respectively. The experimental shapes are obtained from 
the photos using an image analysis program. For vertical jet the only differences between 
the experimental and theoretical shapes is due to bending regions near the nozzle and 
the belt in the experiments. Therefore, we do not discuss vertical jet here. 

In Figure [15] we compare the jet shapes for the concave flow . The experimental shape 
bends steeper down at the nozzle than the one from the model, and is more curved. 
This results in a smaller £ cnc i from the experiment than predicted theoretically. The 
bending region near the belt in the experimental jet does not contribute significantly to 
the difference in x cn( j. 



20 A. Hlod, A.C.T. Aarts, A.A.F. van de Ven and M.A. Peletier 

The comparison of the jet shapes for the convex flow is presented in Figure [16] In the 
experimental shape a small bending region appears near the nozzle where the jet bends 
from the nozzle orientation to some preferred orientation. This orientation is closer to 
the vertical than our theory predicts. The middle part of the experimental jet is almost 
straight, clearly less curved than the one obtained form the model. Near the belt the 
experimental jet bends to the horizontal more rapidly then the one from the model. The 
theoretical shape is curved more uniformly than the experimental one. All this results in 
a larger x on( j predicted theoretically than observed experimentally. 

6.3. Discussion about differences and similarities 

From comparing the jet shapes, including x cn d, for the two flow regimes concave and 
convex, we observe significant differences in the jet shapes obtained theoretically and ex- 
perimentally. The common difference for these flow regimes is the difference in curvature 
of the shapes. Also differences due to bending near the nozzle and the belt are found. 
The steep bending downwards near the nozzle of the concave jet in the experiment, and 
the differences in the jet orientation near the nozzle in the convex flow play a significant 
role in the observed shape differences. As a consequence the theoretical predictions of 
^end are larger than the experimental ones. 

We conjecture that the cause of these differences lies in the effects which we did not 
include in our model, such as air resistance, bending stiffness and surface tension. Ex- 
periments for dnozzic = 1 mm, and L = 0.054m are shown in Figure|2] We found that the 
differences mentioned above in x on d are smaller for the thicker jet falling from a smaller 
height L, which makes us believe that air resistance is important. For the thicker jet, we 
do not observe in the experiments a steep bending of the concave jet near the nozzle; see 
Figures 2(a)[ |2(b)| and [5] The effect of the bending at the nozzle can be compensated 



by adjusting the value of a nozz \ c in our model. Bending stiffness is less important for 
thinner jets, and surface tension for larger flow velocities. For Ubelt close to zero, the 
jet is unstable near the belt, but adequate modeling of this instability is still an open 
question. 

To conclude, we state that our model predicts correctly the transitions between the pa- 
rameter regions for the three flow regimes. Also the tendencies in the (partial) monotonic 
behavior of a; cn d as Ubelt increases are predicted well, yielding a satisfactory qualitative 
agreement. However, significant quantitative differences are obtained. 



7. Summary of the three flow regimes 

Using our knowledge about the three flow regimes from model and experiment, we 
describe typical features of each flow regime. In the model the three flow regimes are 
characterized by the sign of the dimensionless variable £. The value of £ represents the 
momentum transfer through a cross-section of the jet and describes the balance between 
the inertia and viscous terms in the conservation of momentum equation (3.51. Flow 
characterization using experimental jet shape features is possible as well. Below, we 
describe each flow regime separately 

Concave flow. In this flow regime £ is positive. This means that the momentum 
transfer due to inertia is larger than that due to viscosity. This is reflected in the concave 
shape of the jet comparable to a ballistic trajectory. The nozzle orientation is important 
for the jet shape. When the nozzle points vertically down the jet shape in this flow regime 
is vertical, no matter the flow regime is concave or vertical. Therefore, in this case the 
characterization of the flow regime using the jet shape does not distinguish between 
vertical and concave jets. 
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Vertical flow. In this flow regime £ changes sign from negative near the nozzle to 
positive near the belt. Hence, the momentum transfer due to viscosity is larger near the 
nozzle and the one due to inertia is so near the belt. The belt and nozzle orientations are 
now irrelevant for the jet shape, which is straight vertical in the experiments (except a 
possible bending region near the nozzle and bending or unstable region near the belt) as 
well as in the model. 

Convex flow. In this flow regime £ is negative, which means that the momentum 
transfer due to viscosity is larger than that due to inertia. Both in the experiments and 
the model the jet shape is convex (disregarding a small bending region near the nozzle 
in the experiment) and the jet touches the belt tangentially . 

Summarizing, we conclude that the flow regimes can be characterized by the sign of 
the momentum transfer through the cross-section of the jet or by the convexity of the 
jet shape. However, for a n0 zzic = 7r /2 the concave jet shape is vertical, which makes it 
then impossible to distinguish between the concave and vertical flow regimes. Some more 
shape features such as the tangency condition at the belt for the convex flow, and the 
relevance of the nozzle orientation for the concave flow can be used to distinguish these 
flow regimes. 

8. Conclusions 

In this paper we have studied experimentally and theoretically the problem of the fall 
of a viscous jet onto a moving belt. Three flow regimes of the jet are distinguished and 
characterized by the convexity of the jet shape, i.e. concave, vertical, and convex. 

We have modeled the jet using a thin-jet approximation including the effects of iner- 
tia, viscous tension and gravity. The model consists of the stationary conservation laws 
for mass and momentum. A change of the independent variable is made to allow for 
a transformation of the model equations into an algebraic equation. The partitioning 
of the parameter space between the three flow regimes is evaluated in terms of three 
dimensionless numbers. 

The model shows that the sign of the momentum transfer through a cross-section 
of the jet determines the corresponding flow regime. For each flow regime the correct 
boundary condition for the jet orientation is derived by looking at the characteristics of 
the dynamic conservation of momentum equation. These boundary conditions for the jet 
orientation are: 

(a) the nozzle orientation for the concave jet, 

(6) no boundary condition for the vertical jet, 

(c) the tangency of the jet at the belt for the convex jet. 
The missing boundary condition for the vertical jet is replaced by the constraint that at 
the point where the momentum transfer equals zero the jet is aligned with the vertical 
direction of gravity. 

It is shown that a continuous transition between the concave and the convex jets is 
only possible via the vertical one. Also the way how the dimcnsionfull parameters should 
be changed in order to leave the convex or concave jet region is indicated. 

Comparison of the relations between the horizontal position of the touchdown point 
x cnd and the belt velocity Wbcit, obtained from experiments and from the model, shows 
that: 

(a) The model and experiments show similar monotonic behavior of x cnc i as Ubeit is 
changed. 

(b) The parameter regions in the (wbcit , £)-plane for the three flow regimes predicted 
by the model agree with the experimental data. 
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(c) Quantitatively the relations between x cn( j and v^eit show a significant mismatch 
(experiments give smaller x e nd) due to differences in the shapes of calculated and exper- 
imentally observed jets. 

As a final conclusion, we state that the model, which includes viscous tension and in- 
ertia, but disregards air resistance, bending stiffness, and surface tension, describes in 
qualitative sense the fall of a jet of a Newtonian fluid under gravity. 
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